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L INTRODUCTION 

Pulsed inductive plasma accelerators are electrodeless space propulsion devices where a capaci- 
tor is charged to an initial voltage and then discharged through a coil as a high-current pulse that 
inductively couples energy into the propellant. The field produced by this pulse ionizes the propel- 
lant, producing a plasma near the face of the coil. Once a plasma is formed if can be accelerated 
and expelled at a high exhaust velocity by the Lorentz force arising from the interaction of an 
induced plasma current and the magnetic field. While there are many coil geometries that can be 
employed to inductively accelerate a plasma, in this paper the discussion is limit to planar geome- 
tries where the coil take the shape of a flat spiral. A recent review of the developmental history 
of planar- geometry pulsed inductive thrusters can be found in Ref. [1]. Two concepts that have 
employed this geometry are the Pulsed Inductive Thruster (PIT) [2, 3] and the Faraday Accelerator 
with Radio-frequency Assisted Discharge (FARAD) [4]. 

There exists a 1-D pulsed inductive acceleration model that employs a set of circuit equations 
coupled to a one-dimensional momentum equation. The model was originally developed and 
used by Lovberg and Dailey[2, 3] and has since been nondimensionalized and used by Polzin et 
al. [5, 6] to define a set of scaling parameters and gain general insight into their effect on thruster 
performance. The circuit presented in Fig. 1 provides a description of the electrical coupling 
between the current flowing in the thruster Ii and the plasma current / 2 . The acceleration model is 
limited, however, because it lacks a time-accurate description of the energy in the plasma, relying 
instead on an assumed, constant electron temperature T e . While this does permit the computation 
of a plasma resistivity, that resistivity value cannot vary in time. Additionally, various energy sink 
descriptions cannot be included in the model without an equation that accurately describes the 
energy in the system. 

In this paper we revise the pulsed inductive acceleration model to calculate the energy in the 
plasma as a function of time. This will permit the plasma temperature and commensurate pressure 
to vary throughout the duration of the current pulse. It will also provide a path for the inclu- 
sion of more complex energy dissipation models (plasma radiation, two-temperature modeling, 
etc.) without requiring computations through much more complicated two- and three-dimensional 
magnetohydrodynamic simulations. The energy equation is used to compute the deposition of 
power into various modes. One example of these calculations is presented in Fig. 2. The model 


A) 




FIG. 1: A) General lumped element circuit model of a pulsed inductive accelerator. B) Equivalent electrical 
circuit model. (After Ref. [2]) 



FIG. 2: Calculated time-history of the total plasma energy and each component of energy. 


is compared to the case that assumes T e is constant to show the benefits of including the energy 
equation in the calculation. 
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A model for pulsed inductive plasma acceleration containing an energy equation to 
account for the various sources and sinks in such devices is presented. The model consists 
of a set of circuit equations coupled to an equation of motion and energy equation for the 
plasma. The latter two equations are obtained for the plasma current sheet by treating it as 
a one-element finite volume, integrating the equations over that volume, and then matching 
known terms or quantities already calculated in the model to the resulting current sheet- 
averaged terms in the equations. Calculations showing the time-evolution of the various 
sources and sinks in the system are presented to demonstrate the efficacy of the model, with 
two separate resistivity models employed to show an example of how the plasma transport 
properties can affect the calculation. While neither resistivity model is fully accurate, the 
demonstration shows that it is possible within this modeling framework to time-accurately 
update various plasma parameters. 

I. Introduction 

P ulsed inductive plasma accelerators are electrodeless space propulsion devices where a capacitor is 
charged to an initial voltage and then discharged through a coil as a high-current pulse that inductively 
couples energy into the propellant. The field produced by this pulse ionizes the propellant, producing a 
plasma near the face of the coil. Once a plasma is formed if can be accelerated and expelled at a high 
exhaust velocity by the Lorentz force arising from the interaction of an induced plasma current and the 
magnetic field. While there are many coil geometries that can be employed to inductively accelerate a 
plasma, in this paper the discussion is limit to planar geometries where the coil take the shape of a flat 
spiral. A recent review of the developmental history of planar-geometry pulsed inductive thrusters can be 
found in Ref. [1]. Two concepts that have employed this geometry are the Pulsed Inductive Thruster (PIT) 2,3 
and the Faraday Accelerator with Radio- frequency Assisted Discharge (FARAD). 4 

There exists a 1-D pulsed inductive acceleration model that employs a set of circuit equations coupled 
to a one-dimensional momentum equation. The model was originally developed and used by Lovberg and 
Dailey 2,3 and has since been nondimensionalized and used by Polzin et a/. 5,6 to define a set of scaling 
parameters and gain general insight into their effect on thruster performance. The circuit presented in Fig. 
1 provides a description of the electrical coupling between the current flowing in the thruster I\ and the 
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Figure 1. A) General lumped element circuit model of a pulsed inductive accelerator. B) Equivalent electrical 
circuit model. (After Ref. [2]) 


plasma current I 2 . The acceleration model is limited, however, because it lacks a time-accurate description 
of the energy in the plasma, relying instead on an assumed, constant electron temperature T e . While this 
does permit the computation of a plasma resistivity, that resistivity value cannot vary in time. Addition- 
ally, various energy sink descriptions cannot be included in the model without an equation that accurately 
describes the energy in the system. 

In this paper we revise the pulsed inductive acceleration model to calculate the energy in the plasma as a 
function of time. This will permit the plasma temperature and commensurate pressure to vary throughout 
the duration of the current pulse. It will also provide a path for the inclusion of more complex energy 
dissipation models (plasma radiation, two-temperature modeling, etc.) without requiring computations 
using much more complicated two- and three-dimensional magnetohydro dynamic simulations. 

II. Governing Equations 

We proceed with a presentation of the equations governing the response of a pulsed inductive accelerator, 
discussing in turn the circuit model, plasma acceleration model, energy equation, and plasma model. 

A. Electrical Circuit Model 

The lumped-element circuit model of a pulsed inductive accelerator presented in Fig. 1A is redrawn in Fig. 
IB so as to permit the application of Kirchoff ’s law to the two current loops. The external circuit (left side of 
the figure) possesses capacitance C, external inductance Lo, resistance i? e , and acceleration coil inductance 
Lc- The plasma also has an inductance equal to Lc and a resistance R p . The two circuits are inductively 
coupled through the acceleration coil, which acts as a transformer with mutual inductance M. The value of 
M is a function of the current sheet position z. The time- varying behavior of this circuit is governed by the 
following coupled set of first-order ordinary differential equations: 

dh _ VL C + (MIi + hL c ) ( dM/dt ) - I 2 MR P - hR e L c 

dt ~ Lc(Lh + L c )-M 2 

dl 2 _ M (dh/dt) + h (dM/dt) - I 2 R p 

~dt ~ ~L~ C 5 

dV _ _h 

~dt ~ ~C 

where V is the voltage on the capacitor. 


(1) 

( 2 ) 
(3) 
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Following Lovberg and Dailey, 7 the total inductance of the system can be written as 


-^tot {%) — L/q Lp (z) 


(4) 


where L p is the position- varying portion (called herein the plasma inductance), which is equal to the total 
inductance presented to the circuit by the circuit elements on the right side of the dashed line in Fig. IB. 
Based on experimental measurements, it has been found that the plasma inductance as a function of axial 
current sheet position can be modeled as 


L p (z) = L c ( 1 - exp ( - 




( 5 ) 


where zq is defined as the electromagnetic decoupling length. Another way to think of L p is as the approxi- 
mate Thevenin equivalent plasma inductance to the right of the dashed line in Fig. IB. 

The series and parallel inductances in the system can be approximately combined into one total inductive 
element with the value 

/ M 2 \ 

Atot = Ao + f Lc — • (6) 

By equating the terms in the parentheses with the plasma inductance L p , we find that the mutual inductance 
is given as 


M = Lc exp ( — 

The preceding equation can be differentiated to yield 


2 z 0 


dM Lc ( z \ dz 

dt 2 zq ^ \ 2zq ) dt 


( 7 ) 


(8) 


which governs the time- varying behavior of the mutual inductance and completes the set of circuit equations. 


B. Plasma Acceleration Model 

The current sheet, shown schematically in Fig. 2, is modeled as a plasma slab moving with velocity v z into 
an ambient neutral gas. The mass of propellant accumualted by the current sheet is m, the thickness of the 
current sheet is S a , and the temperature, pressure, and ratio of specific heats inside the plasma are given 
as T cs , p cs , and 7 CS , respectively. Ambient neutral gas conditions in front of the sheet are denoted with the 
subscript A and are given for the density pA , pressure pa , temperature Ta, and ratio of specific heats 7 a- 
The sheet is assumed to entrain all the gas it encounters, leaving vacuum conditions behind the plasma. The 
plasma spatially approximated is a disk possessing an inner radius a and an outer radius 6, corresponding 
to the inner and outer radii of the acceleration coil. 

The equations of motion are generally similar to those given in Ref. [7], but 
they have been rederived from first principles, specifically including the pressure 
gradient term in the momentum equation. The equation for mass accumulation 
can be found by assuming one-dimensional motion and integrating the MHD 
equation for conservation of mass over the sheet volume to obtain 


dm . N 


where pA {z) = < ^ 


(9) 


where the assumed mass distribution in front of the current sheet is triagular. 
The linear mass density po corresponding to this distribution is given as 2 (mbit — 
mo ) / dm where mbit is the total propellant mass in the pulse and mo is the initial 
mass in the current sheet. Assuming a value for the ambient temperature Ta, 
the pressure pa can be found using the ideal gas equation of state 



Pa = 


PA 


ks 


T a 


7r (b 2 — a 2 ) M 

where ks is Boltzmann’s constant and M is the atomic or molecular mass 


(10) 


Figure 2. Schematic show- 
ing the variables used to 
describe conditions within 
and just outside the current 
sheet. 
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The momentum equation for the current sheet is found by assuming one- 
dimensional motion and integrating the MHD momentum equation over the current sheet volume, resulting 
in 


dv z 

dt 


~ LcP 

_ 2z 0 


exp 



- p A (z) v\ - pak ( b 2 


/to (t) 


( 11 ) 


where the final term on the right hand side containing pa arises from the pressure gradient term in the 
momentum equation. 


C. Energy Equation 

In previous work the system of equations was closed 7 by assuming that the plasma temperature T was 
constant. The value of this constant was left as a free parameter and varied to produce agreement between 
the model and experimental data. An even more simplified closure approximation that was used in Ref. [5] 
is to assume R p is a constant. 

To properly model the time- varying energy distribution in the plasma we begin with the MHD equation 
governing the evolution of the energy density in the system, given as 8 


de 

dt 


+ V- 


(e + p) v - B m 


= V- 


-E' x B 

do 


( 12 ) 


where e the energy density defined as 


P 1 2 _^ 2 

6 ~ 7 + ~pV + . 

7-1 2 2/i 0 


The first term on the right hand side in the definition of 5 is the internal energy density of the plasma, the 
second term is the gas kinetic energy density, and the third term is the magnetic field energy density. In Eq. 
(12), the term on the right side is Ohmic heating, the term V • (Bm * v'j is the electromagnetic work, the 

term V • (ev) is the net convection of power into a control volume, and V • (pv) represents the pressure work 
performed against the current sheet per unit time. 

We previously observed that the plasma inductance in Eq. 5 is the approximately equal to the Thevenin- 
equivalent inductance to the right of the dashed line in Fig. IB. If we could find a Thevenin-equivalent 
resistance R' p , then the voltage drop across the plasma (i.e. across the two points intersected by the dashed 
line) is 

V p = R P I + 0 (t) = R' p h + (magnetic flux 0 = L P I) . (13) 

Multiplying by I\ yields the total power into the plasma as: 


r d _ (htl\ 

dt V 2 ) 


If dL p 
2 dt 


+ ifRp 


(14) 


where the first term on the right is the rate of change of the magnetic field energy, the second term is the 
rate of electromagnetic work performed, and the third term is resistive heating with the Thevenin-equivalent 
resistive term l\R! p replaced by the actual resistive power load I^Rp- For completeness, the time evolution 
of L p is governed by the equation 


dLp L c ( z \ dz 

~~r~ — — ex P — • 

dt z o \ zq J dt 


(15) 


To use Eq. (12) in the model, the equation must be integrated over the current sheet volume to obtain 
an equation for the total plasma energy E. Doing this transforms the time-differential as 


de dE 
dt dt 


The electromagnetic power term in Eq. (12) is easily replaced using the electromagnetic work term in Eq. 
(14), resulting in 



I\ dLp 

2 dt 


Leif 

2z 0 


exp 



Vz 
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To integrate the term on the right hand side, we must assume that E' = 77 j and that j is curl free, allowing 
this term to transform into the resistive power in the current sheet, 


V- 


-E' x B 
Mo 


I 2 Rp 


The power density convected into the current sheet is 


V • (ev) 


2 P AV Z + 


Pa 

1 A ~ 1 


7 r 


(6 2 



V z 


d f L P I I \ 

dt \ 2 ) 


where the first term on the right represents the power lost in accelerating entrained gas to the sheet speed, 
the second term is the net internal power in the entrained gas convecting into the current sheet, and the 
third term is the time rate-of-change of the electromagnetic field energy as given in Eq. (14). The remaining 
term from Eq. (12) is the work-per-unit-time performed by the ambiant gas pressure against the face of the 
current sheet, and is given as 

V • (pv) -A p A 7T ( b 2 - a 2 ) v z 

Bringing all the transformations above together yields an energy equation for the system as 


dE 

dt 


LJl 

2 z 0 


exp 



V z + H 77 


d 

dt 



2 P AV Z + 


7 aPa 
1A ~ 1 


( b 2 



Vz 


(16) 


After every timestep, the current sheet pressure p cs can be updated using the relation 

E ~ - 7. 1 5 a n (b 2 - a 2 ) + t mv 2 _|_ hlh. 07) 

where 7 CS is the specific heat ratio for the plasma (assumed equal to 5/3). Once p cs is found, the current 
sheet temperature can be approximately computed using an ideal gas equation of state to obtain 


Pcs 


(1 + Z) n e k B 


(18) 


where k B is Boltzmann’s constant, Z is the ionization fraction (equal to one assuming a singly-ionized gas) 
and n e is the electron number density given (for Z = 1) as (m/M)/{5 a tt ( b 2 — a 2 )) where M is the atomic 
or molecular mass of the ions. The temperature is used in the subsequent section to obtain the plasma 
resistivity and other temperature-dependent plasma properties. 


D. Plasma Model 

Following Ref. [7], we can model the plasma as a 
finite width slab that grows in time. The plasma 
resistivity can be simply given using the Spitzer re- 
sistivity, 

rj(t) = 6 x 10 -4 T“ 3 / 2 ohm-m ( T cs in eV). (19) 

An improved (and more realistic) way to find the 
resistivity is to calculate it directly using collision 
cross sections to obtain a more accurate function of 
p as a function of T cs . In Fig. 3 are plotted a Spitzer 
resistivity and an improved resistivity function that 
to an extent captures the non-monotonic behavior 
of an argon gas. We emphasize that the improved 
resistivity curve depicted in the figure is NOT truly 
accurate and is only meant to serve as an example 
of how more realistic resistivity modeling can be in- 
corporated into the system, demonstrating how it 
affects the resulting calculations. 



temperature (eV) 

Figure 3. Comparison of the Spitzer resistivity and an 
improved resistivity model (from Ref. [8] for an argon 
gas. 
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The current sheet thickness S a has an initial, finite value S s at time t — 0 and is then assumed to grow 
on the timescale of resistive diffusion, yielding a time evolution in current sheet thickness as 


Sa(t) = ( (>s + —t 

Mo 


1/2 


( 20 ) 


where po is the magnetic permeability. For a plasma ring having an outer radius 6, an inner radius a, and a 
thickness S a (t), the total plasma resistance is 


Rp(t ) 


7 vrj(t) ( b + a) 


S a (t) ( b-a ) ' 

III. Plasma Simulations 


(21) 


The complete differential equation set given in the previous section consists of Eqs. (1), (2), (3), (8), (9), 
(11), (15), and (16). These differential equations, together with the supporting equations and models for 
temperature, pressure, and resistivity, can be solved for a set of initial conditions. In the following section 
solutions to the equation set are presented for various plasma conditions to demonstrate how the energy is 
partitioned into different modes by the new energy equation. 


improved resistivity model 


Spitzer resistivity 




C) 


3 4 

time (ps) 


D) 


3 4 

time (ps) 


Figure 4. Results of calculations with the acceleration model showing the energy distribution in a pulsed 
inductive system. The results on the left were generated using the improved resistivity function in Fig. 3 while 
the results on the right correspond to the use of a Spitzer resistivity. The mass bits for each case are given. 


In general, the conditions for the simulations are as follows. The electrial properties employed are similar 
to the PIT MkVA. 3 The values of Lq and Lc are 80 nH and 660 nH, respectively, while the value of zo is 
7.5 cm. The capacitance of the bank is 9 /iF and it is charged to 30 kV per pulse. The external resistance 
of the circuit R e is 5 m£2, and the specific heat ratios are all 5/3. The inner and outer radii of the plasma, 
a and 5, are 20 and 50 cm, respectively, and the ambiant gas is assumed to have a temperature of 300 K. 
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The time histories of various calculated energy sinks in a pulsed inductive thruster are presented in Fig. 
4, with results obtained using the improved resistivity model presented in 4A,C and those obtained using a 
Spitzer resistivity in 4B,D. In all cases, the calculation is terminated at the end of the first half-cycle of the 
discharge when the current I\ passes through zero. 

In every graph, the total plasma energy is given as the bold solid line. The total energy is split according 
to Eq. (17) into three components; the internal plasma energy, the electromagnetic field energy, and the 
plasma kinetic energy. The evolution of all three components are given in each graph, summing to the total 
plasma energy. 

We observe that the first half-cycle period is longer for the lower mass bit cases. The internal plasma 
energy is greater for the higher mass bit cases, implying significantly increased heating. In addition, the 
overall energy into the system is greatly reduced in the higher mass bit case owing to poor matching between 
the external driving circuit and the natural acceleration timescale for the propellant (poorly matched dynamic 
impedance 5 ). 

The energy distribution time-histories for the 
two resistivity models are similar, with slightly in- 
creased internal plasma energy at the expense of 
plasma kinetic energy in the Spitzer resistivity case. 

However, from Fig. 5 we find that this difference 
translates into a large difference in the plasma tem- 
perature. The ‘knee’ feature observed in the im- 
proved resistivity model appears to regulate and 
limit the maximum plasma temperature. At lower 
temperatures, resistive heating tends to increase the 
temperature, but as the resistivity drops precipi- 
tously, the plasma heats very little. In the Spitzer 
model, there is no ‘knee’ or inflection point, so 
the resistive heating rate remains high for a much 
greater temperature range until the gas decouples 
from the coil. The initial plasma temperature for 
all these cases was 4.4 eV. However, simulations us- 
ing the improved resistivity model at a mass bit of 
1.2 mg/pulse revealed that after a short transient all 
temperature time histories collapsed to the same curve for an initial temperature range from 0.95 eV to 9.5 
eV. 



Figure 5. Time-evolution of the plasma temperature for 
simulations using the Spitzer and improved resistivity 
models from Fig. 3. The mass bits for each case are 
given. 


IV. Future Work 

It cannot be emphasized enough that while the framework for the energy equation developed in this 
paper is sound, the computations are not yet expected to quantitatively match experimental data. While 
the resistivity model is an improvement over the Spitzer model, it is still not an exact model of the plasma 
resistivity and was used herein to demonstrate that there is an effect on the overall energy distribution in 
the system when a different resistivity model is employed. 

Future work will involve the calculation of several plasma parameters in a more self-consistent and time- 
varying manner. These parameters include the ionization fraction Z, obtained from a Saha relation; the 
plasma temperature, obtained using a real gas equation of state; and the plasma specific heat ratio 7 cs . 
Fortunately, if these data are tabulated then a table lookup can be developed to update these terms in the 
equation set. The resistivity can be found using collision cross section information, calculating it from the 
cross sections for collisions between electrons and every other heavy species in the plasma (neutral, ion, 
double ion, etc.). The completion of this process will permit the customizing of the acceleration model to a 
given propellant type, which is not possible without first having an energy equation incorporated into the 
system of equations. 
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V. Conclusions 


We have presented presented a pulsed inductive plasma acceleration model that includes an energy 
equation to model the various energy sources and sinks in the system. The MHD equations were integrated 
over the plasma volume by assuming one-dimensional fluid motion and treating the current sheet as a 
one-element finite volume. This method generally recovered the equations of motion presented in previous 
work, with the addition of minor corrective terms. More importantly, it permitted the casting of the MHD 
energy equation into terms that were already calculated in the acceleration model, allowing for the time 
evolution of plasma pressure and, through an equation of state, plasma temperature. Results were obtained 
showing the relative sizes of the energy sinks when a Spitzer resistivity and an improved resistivity were 
used, demonstrating that the change in the form of the resistivity as a function of temperature did have an 
effect. While the improved resistivity model is not completely accurate for the plasma, it does show how 
parameters such as the resistivity, ionization fraction, and specific heat ratio, and a real gas equation of state 
can be implemented into future modeling efforts. 
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Circuit Model 



• Pulsed Inductive Thruster circuit coupled to the plasma 


Re Lq Switch 



• Pulsed Inductive Thruster equivalent circuit 



2 


Circuit Equations 



• Current, Voltage equations 

dli _ V L c + ( Mh + I 2 Lc) ( dM/dt ) - I 2 MR p - hR e L c 
dt ~ L c {L 0 + L C ) ~ M 2 

dh _ M (dh / dt) + I\ (d M/dt) - I 2 R P 
dt Lq 

dV _ h 
~dt ~ ~~C 

• System Inductance 

L>tot (^) L 0 + Lp [z'j 

• Mutual Inductance 

„ = iomp (__!_) 
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vacuum 


Plasma Acceleration Model 



current 

sheet 


v z - 

rn., Tcs 
Pcs-) ^cs 


> 


neutral 

gas 

PAi Pa 

Ta, 7 A 


< 


> 


dv z 

dt 


• Continuity Equation 


dm , , 

ip =pa{z)v * 


where Pa(z) ~ 


po (1 Z j Sm) 

0 


Conservation of Momentum Equation 


~7r— exp (~— ) - Pa (z)vl - paw (b 2 - a 2 ) 

ZZq \ Zq 


/m (t) 


Pa 
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Energy Equation Formulation 


Energy density evolution equation 

(e + p) v - Bm ■ v 


de 

h V 

dt 


Energy density 


£ = 


P 


1 


= V 




7 — 1 2 2^o 


-E' x B 
Mo 


Power into plasma from external circuit 

V p = Ftp I + 4>{t) = R'ph + L P ~L + Ii dL 


V 


dt 


dt 


(magnetic flux (f> = L p I) 


P dt 


d ( L v l‘l \ l‘l dL p ( r2 


+ inf + I ! R ” 


Time-evolution of the plasma inductance 

dL p L c f z \ dz 

dt zo eXP I zq) dt 
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Energy Equation Formulation 


de 

dt 


+ V 


(e + p) v - Bm * v 


= V* 


-E' x B 
Mo 


Integrating over the current sheet volume and replacing terms with available quantities 


Time-rate of change of energy 

de dE 

-A 


dt 


dt 


V* 


Resistive heating 
-E' x B 


Mo 


— > I 2 R P 


Electromagnetic work 

-viB U A^-i iLr - LA 


2 dt 


2 zcT eX ^ V 20 


V • (ev) — > 


Gas entrainment 

1 2 

t;Pav z 


Magnetic field energy change 


■ Pa ( j 2 2 \ 

H 7r o — a ) 

1A - 1 V , 

Internal energy flux 


V* 


d_ Lpl{ 
dt \ 2 


Pressure work on control volume 

.2 „2 


V • (pv) ->P47T (6 2 - a 2 ) 
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Energy Equation 



• New energy equation 


dE_ 

dt 


uA 

2^o 


exp 



Vz + “I 77 


d_ 

dt 


T T 2 

2 


2 P* v l + 


7A - 1 



^2 


Total energy in system 

E = - - S a ?r ( b 2 - a 2 ) + ^mt; z 2 + 

7 — 1 2 2 


Pcs c /t 2 „2\ i 1_„.2 i Lpl ^ 


Temperature (ideal gas equation of state) 

Pcs 


T — 

A cs — 


(1 + 2 ^) Tl e kB 
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Plasma Models 



• Resistivity (Spitzer) 

rj(t) = 6 x 10 -4 T“ 3 / 2 ohm- in (T cs in eV) 

• Resistivity (improved model) 

• Caution: not truly accurate, for 
illustrative purposes 




Current sheet growth 

$a(t) = (ds + 

V Mo ) 


1/2 



temperature (eV) 


• Plasma resistance 



7rr/(f) (6 + a) 
(b ~ a) 
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energy (J) > energy (J) 


Plasma Simulations 



c) 


improved resistivity model Spitzer resistivity 




Poorer dynamic 
impedance matching 


More energy into 
internal plasma heating 
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plasma temperature (eV) 


Plasma Simulations 




Much higher temperature for Spitzer 
resistivity model 

Improved resistivity model self-regulates the 
temperature about the ‘knee’ in the resistivity 
vs. temperature graph 

• Colder plasma —*■ higher r| — ► more 
heating 

• Hotter plasma —*■ much lower ri —*■ less 
heating 

Improved resistivity model temperatures 
stabilize to the same value for a wide range 
of initial temperatures (at least 0.95-9.5 eV) 
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Future Work 



• Fully self-consistent resistivity model calculated using collisional cross 
sections for electon impact with multiple heavy species including: neutral, 
ion, double ion, etc. 

• Coupling to models for ionization fraction (Saha), real gas equation of 
state, non-ideal specific heat ratio 

• Permits customization of the acceleration model to a specific propellant 
type 
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